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Abstract 

In the present work, the notion of Super Fractal Interpolation Function (SFIF) is 
introduced for finer simulation of the objects of the nature or outcomes of scientific 
experiments that reveal one or more structures embedded in to another. In the con- 
struction of SFIF, an IFS is chosen from a pool of several IFS at each level of iteration 
leading to implementation of the desired randomness and variability in fractal interpo- 
lation of the given data. Further, an expository description of our investigations on the 
integral, the smoothness and determination of conditions for existence of derivatives of 
a SFIF is given in the present work. 
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1 Introduction 



Barnsley pQ introduced Fractal Interpolation Function (FIF) using the theory of Iterated 
Function System (IFS). Since then, a growing numbers of papers have been published show- 
ing relation between fractals and wavelets jH [10], fractal functions and Kiesswetter-like 
functions [2] and on fractal dimension [131 E]. Later, Barnsley et. al. extended the idea 
of FIF to produce more flexible interpolation functions called Hidden-variable FIF (HFIF) 
which were generally non-self affine. Dalla [7] found bounds on fractal dimension for the 
graphs of non- affine FIFs. In 1989, Barnsley and Harrington [4] constructed an IFS to show 
that a FIF can be indefinitely integrated, giving rise to a hierarchy of smoother functions 
and developed results on differentiability of a FIF. 

Fractal Interpolation Function, constructed as attractor of a single Iterated Function System 
(IFS) by virtue of self-similarity alone, is not rich enough to describe an object found in 
nature or output of a certain scientific experiment. The objects of nature generally reveal 
one or more structures embedded in to another. Similarly, the outcomes of several scientific 
experiments exhibit randomness and variation at various stages. Therefore, more than one 
IFSs are needed to model such objects. Barnsley [HOE] introduced the class of super fractal 
sets constructed by using multiple IFSs to simulate such objects. Massopust p2] constructed 
super fractal functions and V-variable fractal functions by joining pieces of fractal functions 
which are attractor of finite family of IFss. However, for a data set arising from nature 
or a scientific experiment, a solution of fractal interpolation problem based on several IFS 
has not been investigated so far. To fill this gap, the notion of Super Fractal Interpolation 
Function (SFIF) is introduced in the present work. The construction of SFIF requires the 
use of more than one IFS wherein, at each level of iteration, an IFS can be chosen from 
a pool of several IFS. This approach is likely to ensure desired randomness and variability 
needed to facilitate better geometrical modeling of objects found in nature and results of 
certain scientific experiments. The construction of SFIF is followed in the present paper by 
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the investigations of its smoothness, its integral and determination of conditions for existence 
of its derivatives. 

The organization of the present paper is as follows: In Section |2j for a given finite set 
of data, the method of construction of a Super Fractal Interpolation Function (SFIF) is 
developed. At each level of iteration, an IFS is chosen from a pool of IFS in our construction 
of SFIF. For a sample interpolation data, a computational model of SFIF, illustrating the 
construction method given in Section [2} is presented in Section [3j The fractal dimension and 
average fractal distance are computed for various SFIFs constructed in this section. Finally, 
in Section HI it is found that for a SFIF passing through a given interpolation data, its 
integral is also a SFIF, albeit for a different interpolation data. An expository description 
of smoothness of a SFIF and conditions for existence of derivatives of a SFIF is also given 
in this section. 



2 Construction of SFIF 

In this section, the notion of Super Fractal Interpolation Function (SFIF) is introduced via 
its construction based on more than one IFS. 

Let Sq = {(xi,yi) G M 2 : i = 0, . . . N} be the set of given interpolation data. For k = 
1, . . . , M, M > 1 and n = 1, . . . , N : let the functions u Ui k :/xt- )■ / x M be defined by 

u n ,k{x, y) = (L n (x), G n>k (x, y)) for all 0, y) G M 2 (2.1) 
where, the contractive homeomorphisms L n : I — > I n are given by 

T / \ .7 \X n X n —i)X -\- yX]yX n —i XoX n ) . . 

L n (x) = a n x + b n = r (2.2) 

{x N - x ) 
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and the functions G n ^ :/xl4l defined by 

G n ,k(x, y) = e n:k x + j n>k y + f n>k (2.3) 

satisfy the join-up conditions 

G n ,k{xo,Vo) = Vn-i and G ntk (x N , y N ) = y n . (2.4) 

Here, are free parameters chosen such that |7 ni fc| < 1 and ■jn^ ^ J n ,i f° r k ^ I. By (12. 4p . 
it is observed that u n> k are continuous functions. The Super Iterated Function System 
(SIFS) that is needed to construct SFIF corresponding to the set of given interpolation data 
So = {(xi, yi) G M 2 : i = 0, 1 . . . , N} is now defined as the pool of IFS 

{{R 2 ;a; n , fc :n=l,...,iV}, fc = l,...,Af} (2.5) 

where, the functions u ritk are given by (12. ip . 

To introduce a SFIF associated with SIFS (ESI), let {W k : H(R 2 ) "H(K 2 ), k = 1, . . . , M}, 

N 

be a collection of continuous functions defined by W k (G) = |J uj n ^{G) where, u njk (G) = 

n=l 

u n)k (x,y) for all (x, y) G G. Since, h(W k (A), W k (B)) < max j Utk h(A,B), where h is 

l<n<N ' 

Hausdorff metric on "H(IR 2 ), {H{R 2 ); W\, . . . , Wm} is a hyperbolic IFS. Hence, by Banach 
fixed point theorem, there exists an attractor A G "H("H(IL 2 )). 

Let A be the code space on M natural numbers 1, 2, . . . , M . For a = o\02 ■ ■ ■ cr k . . . G A, 
define the function : A — > "H(IR 2 ) by 

0(a) = lim W Uh o W Gh _ x 0...0 W ai (G), G G H(R 2 ). (2.6) 

It is shown that [2] <fi(cr) exists, belongs to A and is independent of G G "H(M 2 ). Also, the 
function is onto and continuous [2]. In the construction of SFIF, for a a = <J\Oi ... G A, 
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let the action of SIFS (I2.5P at the iteration level j be defined by Sj = W aj (Sj-i), where So 
is the set of given interpolation data. It is easily seen that the set, 

G a = 0(a) = lim W ah o . . . o W ai (S ) = lim S k (2.7) 

fc— >oo k— >oo 

is the attractor of SIFS ( 12.5ft for a fixed a G A. The following theorem shows that G a is the 
graph of a continuous function g a . 

Theorem 2.1 Let G a be the attractor of SIFS (12.51) for a = oyOi . . . a fe . . . G A. Then, G a 
is graph of a continuous function g a : / — > R such that g a {x n ) = Dn for all n = 0, . . . , N. 

Proof Let go be a function whose graph is So- Then, the set S^, k > 1, is graph of the 
function g ak , where g ak (x) = G ik ^ k {L^{x),g ak _ 1 {L^{x))). It is easily seen that g ak (x) = 
Gi k ,a k ^V), £4-1,^-1 G iuai (Lr 1 o...oL; h 1 (x),go(L; i 1 o...oLr 1 (x))) . . . jj . There- 
fore, it follows by (12. 7p that the set G a is graph of the function g a = lim g Uk . 

k—^oo 

For proving the continuity of the function g a , consider r^r^ . . . r* . . . G A where r* ^ 1 for 
some j G N and r* = 1 for i G N and i ^ j. We first show that G T * is graph of a continuous 
function g T *. If not, then G T * = 0(t*) is graph of a function g T * that is not continuous so 
that there exist a 5± > such that whenever Xi,X2 G / and |xi — 22I < Si, 

\g T *{x 1 )-g T *(x 2 )\>e. (2.8) 

It is known that [T], for r = 1 G A, G T = 0(r), with defined by (I2.6p . is graph of a 
continuous function gv : / — > M such that g T (x n ) = y n , n = 0, 1, . . . , N. Consequently, there 
exists a 62 > such that \x\ — x<i\ < 82 implies \g T (xi) — g r (x 2 )| < §. Also, since is a 
continuous map, there exists 5 3 > such that, for r and r* satisfying d c (r, r*) = < 83, 

max|g r (x) — g T *(x)\ < f. Thus, for 5 = min(5 1 , <5 2 , 5 3 ) and Xi,x 2 satisfying \xi — x 2 \ < 8, 
\gr*{xi)-g T *{x2)\ < |fl , r*(^i)-5 , r(^i)|+|fl , r(^i)-fl , r(^2)|+|5'r(^2)-fl , r*(^2)| < e, a contradiction 
to (12.81) . Hence, G T * is graph of continuous function g T * . 
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Now, consider the sequence o~ n = o~\, n <32,n with a^ n = o~j for j < n and o~^ n = 1 for 

j > n. It is easily seen that as n tends to infinity, a n tends to a with respect to the metric 
d c . Using the arguments of previous paragraph inductively, it follows that G an = 4>(a n ) is 
graph of a continuous function g an defined on I. Let G a = 0(cr) be graph of a function 
g a . By continuity of <p, G Un tends to G a with respect to Hausdorff metric h as n — > oo, 
which implies that g„ n tends to g a with respect to Maximum metric as n — > oo. Hence, 
there exist an e > such that max\g crn (x) — g a (x)\ < |. Since g an is continuous on /, 
there exits a 5 > such that \x — y\ < 5 implies \g Un (x) — g an (y)\ < §■ Therefore, 

\g*(x) - g<r(y)\ < \g (x) - g an {x) \ + \g an {x) - g<r„(y)\ + \9a n (y) - g*(y)\ < e for \x - y\ < 5 

implying that the function g a is continuous on /. This establishes that the attractor G a of 
SIFS ( 12. 5 p is the graph of continuous function g a . ■ 

Theorem 12.11 is instrumental in defining a SFIF associated with SIFS ( 12. 5ft as follows: 

Definition 2.1 The Super Fractal Interpolation Function (SFIF) for the given 
interpolation data {(xj, yi) : i = 0, 1, . . . , N} is defined as the continuous function g a whose 
graph G a is the attractor of SIFS (12. 5 j) . 

Remark 2.1 Consider the family of continuous functions Q = {/:/—>■ R such that 
f is continuous, f(xo) = yo andf^x^) = Z/tv} with metric dg(f, g) = max \f(x)— g{x)\. Since 
Q is a complete metric space, it is easily seen that, for a fixed a G A, Read-Bajraktarevic 
operator T : A x Q — > Q defined as 

T(a, g)(x) = fc lim Ig^ (l^( x ), G^^ (Lr 1 ^ o Lr\x), G ik _ 2>ak _ 2 (.,... 

G il>ai (Lri _ _ _ Q L r k \ x ) jg (L^ o . . . o L r\ x ))) -..)))}» ( 2 - 9 ) 

is a contraction map on Q and so it has a unique fixed point in Q . It is observed that, SFIF 
g a satisfies g a = T(a,g a ). 
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Remark 2.2 The notion of SFIF can further be generalized by introducing a fixed parameter 
k (0 < k < 1) in the join-up conditions ( 12.4ft as follows: 

G n , k (xo, kx + (1 - K)y ) = Kx n -i + (1 - K,)y n -i 

G H: k(xN, KX N + (1 - K,)y N ) = KX n + (1 - n)y n . 

The above condition ensures that there exits a unique attractor G CT K G "H(R 2 ) of SIFS ( 12. 5p . 
i?y t/ie arguments similar to those in the proof of Theorem \2.1\ G a>K zs graph of a continuous 
function g aK , called henceforth Parameterized SFIF or n-SFIF. 



3 Computational Model of SFIF 



Our method of construction developed in Section [2] is employed in the present section for gen- 
erating various SFIF for a sample interpolation data So = {(0, 0), (30, 90), (60, 70), (100, 10)}. 
For identifying the corresponding SIFS j{R 2 ;a; nj fc : n = 1,2,3}, k = 1,2|, the maps 
&n,k, k — 1,2 (c.f. ( 12. lft ) are obtained by computing (c.f. Table []]) the values of a^, 
h (c.f. (JZ2D) and e itl , / M ; e ij2 and f i:2 (c.f. (JH) with 7i>1 = 0.4 and 7i , 2 = 0.6 for % = 1, 2, 3. 

In the construction of SFIF for a a = <Ti<t 2 . . . G A, the set Sj = W aj (Sj_i), j = 1,2,..., 
representing the action of SIFS (12.51) at the iteration level j is computed. The SFIF g a for 



cr( b ) = 1 (c.f. Figs. 1(a)- 1(c), blue curve) is constructed by the action of IFS {M 2 ;w ni i,n 



1, . . . ,N} at every level of iteration. Similarly, SFIF g a{g) for cr( 9 ) = 2 (c.f. Figs. 1(a)- 1(c) 
green curve) is constructed by the action of IFS {M 2 ; uj n< 2, n = 1, . . . ,N} at every level of 



iteration. The SFIF g CT(r) for <7( r ) = 112 (c.f. Fig. 1(a) , red curve) is constructed by the action 



of IFS {M 2 ;w n ,i,n 



N} at j th level of iteration if j is not divisible by 3 and by the 



action of IFS {M 2 ; u n> 2, n = 1, . . . , A^} if j is divisible by 3. Likewise, SFIF g a{r) for = 221 



c.f. Fig. 1(b), red curve) is constructed by the action of IFS {R ,ui n i,n 



N} at j 



tii 



level of iteration if j is divisible 3 and otherwise by the action of IFS {R ; u n> 2, n 



N}. 
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Finally, SFIF g a for <j( r ) = 12 (c.f. Fig. 1(c) , red curve) is constructed by the action of IFS 
{IR 2 ; w n ,i> n = 1, . . . , N} at j th level of iteration if j is not divisible by 2 and by the action of 
1, . . . , N} if j is divisible by 2. 



IFS {R 2 ;u n , 2 ,n 





i=l 


i=2 


i=3 


a 


0.3 


0.3 


0.4 


b 





30 


60 




0.86 


-0.24 


-0.64 


fi,l 





90 


70 


ei,2 


0.84 


-0.26 


-0.66 


fi,2 





90 


70 



Table 1: Computed Values of a,, 6j, e^i, f^i, e ij2 , fi,2, i = 1, 2, 3, for sample data S 




(a) cr (r) = 112 (b) cr (r ) = 221 (c) cr (r) = 12 



Blue Curve ( — ): SFIFg,j (b) for er( b ) = 1 in(a), (&) and (c), 
Green Curve ( — ): SFIFg a{g) for CT( 9 ) = 2 in(a), (&) and (c), 

Red Curve ( — ): SFIFg a(r) for cr( r ) = 112 in(a), for cr( r ) = 221 in (b) and for er( r ) = 12 in (c) 



Figure 1: SFIFs for = 1, = 2 and different choices of <7( r ) 



The SFIFs g a for o"( 6 ) = 1 and g a for <j( 9 ) = 2 are in fact FIFs (c.f. Figs. 1(a)- 1(c), blue 



and green curves ), since these are constructed with a single element of SIFS (12. 5p . 

Heuristically, in terms of their fractal dimension [T] , the graphs of SFIF g a{r) appear to fill 
more space in M 2 than the graph of FIF g a{h) and less space in R 2 than the graph of FIF g a(g) . 
In fact, the fractal dimension of graphs of FIF g a{b) and g a{g) are computed as 1.3069 and 



la 



1.5199 respectively whereas the fractal dimension of SFIF g a{r) with <j( r ) = 112 (c.f. Fig. 
red curve) is 1.3632, the fractal dimension of SFIF g a[r) with = 221 (c.f. Fig. 1(b), red 



curve) is 1.4572 and the fractal dimension of SFIF g a with = 12 (c.f. Fig. 1(c), red 
curve) is 1.4182. 

Further, for FIFs g a(b) and g a , the average fractal distance defined as 

d>F{f,g) — (bZ^y ( / \ f( x ) — g(x)\ 2 dx J for the functions / and g, continuous on a closed 
interval [a, b], is dF(ga (b) , 9a (g) ) — 0.297. It is observed that (i) for SFIF g a ,. with <7( r ) = 112, 
dp{9u (b y 9a [r) ) = 0.022 while dp(g a , g) , g a{r) ) = 0.276. So, if the data generating function is at 
one third average fractal distance from FIF g a(b) , then SFIF g a[r) is a better approximation 



of the data generating function, since g a{r) is closer to g a(b) than g a , g) (c.f. Fig. 1(a)) i.e. 
dF(ga ib) ,g<T ir) ) < d F (g a(g) ,g a{r) ). (ii) For SFIF g a(r) with a {r) = 221, d F (g am , g a{r) ) = 0.228 
while dp{ga (g) , 9a {r) ) = 0.071. So, if the data generating function is at one third average fractal 
distance from FIF g a(g) , then SFIF g a is a better approximation of such data generating 



function, since g a , . is closer to g a . . than g a (c.f. Fig. 1(b) ) and (iii) for g a . with a^) = 12 



dF(g a(b) , g U(r) ) = 0.138 and dF(ga {g) , 9a (r) ) = 0.162. So, if the data generating function lies in 



the middle of FIFs g a , b) and g a{g) ^ then SFIF g a , (c.f. Fig. 1(c)) is a better approximation 
of such data generating function. 



4 Integral and Derivative of SFIF 

In this section, for a SFIF passing through a given interpolation data, its integral is shown 
to be also a SFIF, albeit for a different interpolation data. Further, in this section, the 
smoothness of SFIF is investigated in terms of its Lipschitz exponent and it is found that, in 
general, a SFIF may not be different iable. This, as a natural follow up, led to determining 
in this section the conditions for existence of derivatives of a SFIF. 

In order to study the integral of a SFIF, a SIFS 

{{M 2 ; u ntk (x,y) = (L n (x),G nik (x,y)):n = l,...,N}, fc = l,...,Af}, (4.1) 
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associated with the data !/,) el 2 :i = 0,..., N} is considered, where L n (x) = a n x + b r , 
are given by (12. 2 \ and the functions G njk (x,y) defined by 



G n>k (x, y) = j n7 kV + q n ,k( x )i n = l,...,N. 



(4.2) 



satisfy the join up conditions given by (12.41) . Here, 7^ are free parameters chosen such that 
|7njt| < 1 and j n;k 7^ j n j for k 7^ I and q n}k (x) are affine functions. Condition (12 .4p ensures 
that there exits a unique attractor G a G ^(IR 2 ) of SIFS (14. ip . By the arguments similar to 



those in the proof of Theorem 12. 11 G a is graph of a continuous function g a . 

The following notations jl] are needed in the sequel for tidy presentation of our results: 



7n,k 
VN,k 



Q"nTfn,k 



N r X N -1 

E a j i [ / <?j,fc(*) d q 



yo + 



jv r X JV 

E 

£=1 3:q 
N 



1_ E a j7j,fc 



^0 + E a i [7j,fe(&v,fc - J/o) + / dt] 



yn,fc 

<?n,fc(z) = yn-i,fc - anln,ki)o + a n J q n , k (t) dt 



(4.3) 



where, ?/o is an arbitrary real number. To determine an interpolation data through which 
the integral of SFIF passes, let the affine functions q nk (x) in (14. 2p satisfy : 



N x N N x N 

E "j I Qj,k E "j I Qj,i 

i =1 XQ 3= l X Q 



N N 

1 - E a iij,k 1 - E a jij,i 



^ 1 for k y£ I, k,l = 1,...,M. 



(4.4) 



For example, for aj = — , j jjk = j k and q jtk = (1 - ^ k ){ejX + j)) for j = 1,...,N, the 
condition (14. 4p is satisfied. Then, y i>k = y^ = y^ for i = 0, . . . , N; k,l = 1, . . . , M and 

The SIFS associated with the data {(xj,^) G M 2 : z = 0, . . . N} is now defined as the pool 
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of IFS 



{R 2 ; u n , k (x,y) = (L n (x),G n 4x,y)): n = 1, . . . , N}, k = 1, . . . , m} (4.5) 



where, the functions 



Gn,k(x, V) = ln,kV + q n ,k( x ) ( 4 - 6 ) 

satisfy the join-up conditions G n ,k(xo,y ) = y n -i and G n ^{x^,y^) = y n . These join- 
up conditions ensure that there exits a unique attractor G„ G "H(K 2 ) of SIFS (14.51) . The 
following theorem shows that the integral of SFIF is also a SFIF albeit for interpolation data 
{{x^y,) el 2 :i = 0,...N}. 

Theorem 4.1 For the interpolation data {(xi,yi) G M 2 : % = 0,...N}, let g a be SFIF 
corresponding to SIFS (14.11) for o G A. Then, the integral 

X 

9a(x) = y + J 9*(t) dt (4.7) 

•I'D 

is SFIF associated with SIFS (14.5p for the interpolation data {(xi, iji) : % = 0, . . . , N}. 
Proof Using ( 14. 7ft and ( 12. 9ft . it is observed that, 



9a(L ih 0...0 L h (x)) = g a (L ik 0...0 L h (x )) + I J]/' ~ ,, J (g a \ 



x) - y 



£ jj) _ 1 o...oL< 1 (s) 



k / k \ 

+ ( n ) % / 



fcp.<r„(*)- (4. 
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Also, by gZD, 



g a (L ik 0...0 Lt^xo)) = yo + J2[ II a H^w ) | 

p=i \j=p+i ' ^ 1=1 

X N 



a i 



ll,a p {VN - Vo) 



J-o 



■t'o 



The above identity and (14. 3[) give 

9a(L ik o . . .oL h (x )) 
k . 

3=1 7 

• n ( n s^,-, 

P=l 3=P+1 



o... 0^(10) 



. (4.9) 



•I'd 



Now, substituting the value of g a (Li k o . . . o (x )) from (HSJ) in p~8|) . it follows that 
g a (L ik 0...0 L h (x)) = G ik)(7k 0...0 L h (x), 



^iitOi f-^ii (^0 ) (*^) fl'c(^'))) • • 



Thus, p CT is SFIF associated with SIFS (14. 5ft . 



Remark 4.1 In case of k-SFIF g CT>K (c.f. Remark \2.2\) , using the lines of proof of Theo- 
rem \4-l\ it follows that the integral of k-SFIF is not a n-SFIF but integral of g atK + £ a , where 
£ CT is defined by 



6x(x) = 7in,«r„ [UL-\x)) - «(1 - L~ 1 (x))] + k(1 - z) 



(4.10) 
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for x G I n , is a k-SFIF for the interpolation data {(£«, yi) G M 2 : i = 0, . . . N}, provided 



N x n N x N 

E a AI - «7j,fc(l - E %[/ - «7j,z(l - t))dt] 



N 



N 



^ 1 holds. 



1 - E «j7j,fc 



1 - E a i7i,i 



i/ere, fc , <^ fc (x), j = 1, 2, . . . , N, in (SI3D are azven by 



(1 - K)?Mr,fc = k(x q - Xjy) + (1 - K)y 



+ 



N f %N X N %N 

E a n\ J Qn,k(t) dt + K f (1 - L n (t)) tZt - « 7 n ,fc / (1 - t) dt 
n=l L x so xo 

TV 

1 - E °n7n,fc 
n=l 



(1 - K)y jjk = k(x -Xj) + (1- n)y Q + ^ a„ 7„ J (1 - k 

n=l L ^ 



0(&v,fc - £o) 



,r A r 



+ Kjtdtj + J {«(1 - L„(t)) + g n , fc (t)} dt 

x icq 



and 



9j,fe(^) = + (1 - «)%-i,fe - aj7j,fc 



KXq + (1 — «)$o + k I (1 — t) dt 



X X 



(4.11) 



•f'O 



■!■() 
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For investigating the smoothness of a SFIF, the following notations are needed: 



A = min{A n> fc : n = 1, 2, . . . , N, k — 1,2, ... M}, where are real numbers 

satisfying < A ni fc < 1 
C\ = max{ 7rf x : n = 1, 2, . . . , TV, k = 1, 2, . . . , M}, where 7 n> fc are real numbers 

I n\ 

satisfying \^ k \ < 1. 

Modulus of continuity of g a (x) as Lo(g a , t) = max max \g a {x + h) — g a {x)\ (4.12) 

\h\<t x 

The smoothness of a SFIF in terms of its Lipscitz exponent is given by the following theorem: 

Theorem 4.2 Let g a be a SFIF corresponding to SIFS (14. ip with q n ^ G Lip X nt k, < \ rh k < 
1 . Then, 

(i) for Ci <1, g a E Lip A 

(ii) for d = l, oo(g a ,t) = 0(\t\ x \og\t\) 
(Hi) for C\ > 1, g a G Lip A, 

where, A < max ^ -^~ J an d C\, A are given by (I4.12p . 

fc=l,'...,M 

Proof The method of proof is similar to that in [9], wherein 7 n is replaced by 7n )(Tn . ■ 

Remark 4.2 It follows from Theorem \4-2\ that g a G Lip A& /or Ci > Ci^ > 1, where 
C l!fc = max{^4 :n=l,2,...,iV} and A^ < max{ig^ : n = 1, 2, . . . , N}. 

Remark 4.3 In case of k-SFIF g aK (c.f. Remark \2.2) . the smoothness result analogous to 
Theorem \4-2\ can be obtained as follows: (i) g aK G Lip A for C\ < 1, (ii) u(g a)K ,t) = 
0(|t| A log \t\) for Ci = 1 and (iii) g a ^ K G Lip A for C 1 >1,\< max{(log7 nifc )/(k>g a n) : n = 
1,2,. ..,N, fc = l,2,...,M}. 
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In general, a SFIF belonging to certain Lipschitz class, need not be differentiable. This, as 
a natural follow up, leads to identification of conditions for the existence of derivative of a 
SFIF in the following proposition: 



Proposition 4.1 For the interpolation data {(x{, yi) G 1R 2 : % = 0, 1, . . . N}, let g a be a SFIF 
corresponding to SIFS (14. ip for a e A. Then, g' a exists and g' a (x) = g a {x) if and only if g a 
is a SFIF associated with SIFS (14. 5 p for the interpolation data {(xi,yi) : i — 0, 1, . . . ,N}, 
provided % k = ctjlj,k and ^(q jjk (x)) = a^q^x) hold. 

X 

Proof If g' a (x) = g a (x), then g a (x) = yo + J Qaif) dt, so that " if " part follows from The- 

x 

orem l4~Tl Conversely, suppose g a is a SFIF associated with SIFS ( 14.51) for the interpolation 
data {(xi, y~i) : i = 0, . . . , N}. Then, 

g a (L ik o . . . o L h (x)) = (l[[%,a^J 9*{x) 

k / k \ 
+ E II %W <W L *p-i ° • • • ° L n(^))- (4-13) 

p=l \j=p+l / 



Since, ^(g jjfc (a;)) = a^q^x) 



dj,k(x) = qj,k( x o) + o-j j qj,k(t)dt = jjj-i - aj>y j:k yo + a, y q jjk (t)dt. (4.14) 
Substituting (14.14ft and 7^ = %7j,fc in (14. 13ft . 

k k k 

g a (L ik 0...0 L h (x)) = (Y[ai.% h(T /)g a (x) + E ( II "■'■,-■) x 

j=i p=i j=p+i 

L ip _ 1 o...°L il (x) 

Vi P -i - a ip ^ p ,a p yo + a ip J q ip ,* p (t)dt . (4.15) 



x 



■!■() 
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For a fixed a G A, it is easily seen that the Read-Bajraktarevic operator T defined by 



{k k k 

( n ">■">" )-' /(/ -' 1 ° • • • ° L 7k( x )) + yi ( n a ^iw 
j=i p=\ j=p+i 



X 



LJ 1 o...oLJ 1 {x) 



X 



J/ip-i - ai p li p ,a p yo + a ip I qi p ,a p (t)dt \ (4.16) 



■I'D 



is a contraction map on £y = {/ : / — > R such that / is continuous, f(xo) = Vo and 
/(xjv) = 2/at}- By (14.151) . the function g CT is a fixed point of T. Also, by Theorem 14. 1[ 

X 

the function h(x) = y + J g a {t) dt is a SFIF associated with SIFS (14. 5p satisfying (14.15)) . 
Consequently, h also is a fixed point of T. Hence, by uniqueness of fixed point of Read- 

X 

Bajraktarevic operator T, g a (x) = yo + / 9a{t) dt which implies that g' a exists and g' a (x) = 

x 

g a (x), since g a being a SFIF corresponding to SIFS (14. ip . is a continuous function. ■ 
For the investigation of n th derivative of SFIF, denote 

G ijk ,j(x, y) = 7i )fc j y + qi, k ,j(x) (4.17) 

where, G itkfi (x,y) = G i)k (x,y), qi, k ,o(x) = Qi,k( x )> 7i,fc,o = Ji,k and Gi^ k> j{x N ,y N>k j) = 
Gi >k> j(xo, yo,k,j), i — 1, ■ ■ ■ , N, k = 1, . . . , M and j — 0, 1, . . . , n. To determine interpolation 
data through which derivatives of SFIF passes, let the affine maps qi }k j(x) in ( 14.171) satisfy : 

N x N N x N 

^2 a p I ip,k,j a p I ip,i,j 

P=l XO P =l X (41g) 



N N 
p=i p=i 

where #o,j>.7 = 0, 1, . . . , ?7- are arbitrary real numbers. For example, for aj — j^, 7i,fc,,7 — 
and g^jfcj = (1 — 7^)^(2;), where q~i,j(x) are polynomials of degree n — j for z = 1, . . . , N, 
the condition (I4.18P is satisfied. Then, y i k j = jjijj = y^ for i = 1, . . . , N, k,l = 1, . . . , M 
and j = 0, 1, . . . , n. 
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The SIFS associated with the interpolation data {(xi, y^j) : i = 0, 1, . . . , iV}, j — 0, 1, . . . , n, 
is now defined as 

{{K 2 ; u ilk>j {x,y) = (L i (x),G i>k j{x,y)):i = l,...,N}, fc = l,...,M}. (4.19) 

It is observed that SIFS ( 14. 19ft reduces to SIFS ( 14.11) if j = 0. The following theorem gives 
the existence of derivatives of a SFIF. 

Theorem 4.3 Let the functions Gi jk) j(x,y) defined in (14. 17ft be such that, for some integer 
n > 0, \^ i>k \ < all, q i>k E C n [x , x N ], i = 1, 2, . . . , N, k = 1, 2, . . . , M and g a be a SFIF 
corresponding to SIFS (14 . 1 9 [) for j = and a G A. Then, for j — 1, 2, . . . , n, g& exists and 
is a SFIF associated with SIFS (14. 19ft for the interpolation data {(xi, yij) : i — 0, 1, . . . , AT}, 
provided j ijkJ = ^ and qi,k,j( x ) = ■ 



Proof The equation Gi^j (x , y j) = y j gives y j = -^rUoj + <?1 '^f°' > which implies 

= 7T v Similarly, G N>k ,j {x N , y NJ ) = y NJ gives y N>j = — f -. By Proposi- 

[ a l ^fl,k) \ a N *yN,k) 

tion 14.11 it now follows that, for j = 1,2, ... ,n, ga\x) is the SFIF associated with SIFS 
{M 2 ; ou iAj (x,y) = (Li(x),G iAj (x,y)) : % = 1, 2, . . . , N], k = 1,2,...,m). ■ 



Remark 4.4 In the case of k-SFIF, Remark \4-1\ and Theorem \4-3\ suggest that the function 



gal = ga,l(x) — $ ^(x), with £ CT given by (14.10)) . is k-SFIF associated with the SIFS 
P =i 

{{M 2 ; u iikJ :i = 1,2,..., N}, k = 1, 2, . . . , m}, j = 0, 1, . . . , n. 
Using (14. 10p and (14 .11)) . it is easily seen that y ,i = ( - 1 _^' ( -^ i _ 7i n — jz^ and 

1 { Nk( x N) K , , <llk( x 0) , 

VNki—n — w \ — -r^- Also, Voi = -, rn r ™«jii,- = : -j-, -. 

fc = l,2,...,M,/orj>l. 
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5 Conclusions 



In the present work, the notion of Super Fractal Interpolation Function (SFIF) is introduced 
for finer simulation of the objects of the nature or outcomes of scientific experiments that 
reveal one or more structures embedded in to another. Since, in the construction of SFIF, 
at each level of iteration, an IFS can be chosen from a pool of several IFS, the desired 
randomness and variability can be implemented in fractal interpolation of the given data. 
Thus, SFIF may be used as a tool for better geometrical modeling of objects found in nature 
and results of certain scientific experiments. Also, an expository description of investigations 
on the integral, the smoothness and determination of conditions for existence of derivatives 
of a SFIF is given in the present work. It is proved that, for a SFIF passing through 
a given interpolation data, its integral is also a SFIF, albeit for a different interpolation 
data. The smoothness of a SFIF is given in terms of its Lipschitz exponent. A SFIF 
g a , for C\ 7^ 1, belongs to a Lipschitz class and, for C\ = 1, oj{g a ,t) = 0(\t\ x log \t\). 
It is seen that the smoothness of SFIF depend on free variables 7 n> fc as well as on the 
smoothness of affine functions q n ,k{x) occurring in its definition. Further, sufficient conditions 
for existence of derivatives of a SFIF are derived in the present paper. Our results on SFIF 
found here are likely to have wide applications in areas like pattern-forming alloy solidification 
in chemistry, blood vessel patterns in biology, signal processing, fragmentation of thin plates 
in engineering, stock markets in finance, wherein significant randomness and variability is 
observed in simulation of various processes. 
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